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Abstract 



We match collapsing inhomogcneous as well as spatially homogeneous 
but anisotropic spacetimes to vacuum static exteriors with a negative 
cosmological constant and planar or hyperbolic symmetry. The collapsing 
interiors include the inhomogeneous solutions of Szekeres and of Barnes, 
which in turn include the Lemaitre-Tolman and the McVittie solutions. 
The collapse can result in toroidal or higher genus asymptotically AdS 
black holes. 



1 Introduction 

The process of gravitational collapse to a black hole in general relativity can 
be modelled by the matching of two appropriate exact solutions. For a zero 
cosmological constant A there are several well-known examples in spher- 
ical symmetry, one of the simplest being the Oppenheimer-Snyder model 
which results from the matching of a Friedmann-Lemaitre- Robertson- Walker 
(FLRW) dust metric to a Schwarzschild exterior. For A / spherical cases 
the results may be less well-known: the matching of A-FLRW with Kottler 
has been done by Balbinot et al. [1], Nakao [19] and Markovic & Shapiro 
[18], and its A-Lemaitre-Tolman counterpart by Lake [12]. The matching of 
the dust A-Szekeres solution with Kottler was recently studied by Debnath 
et al. [7] (following Bonnor [4] for A = 0), while the matching of a collapsing 
fluid with tangential pressure and A ^ to Kottler has been investigated by 
Madhav et. al. [17]. The above examples give rise to spherical black hole 
solutions with A^O. 

The collapse to non-spherical black holes is less widely studied. Smith 
<fc Mann [20] have shown that one can match a collapsing k = — 1 FLRW 
spacetime to an asymptotically anti-de Sitter (AdS) exterior, as a model of 
gravitational collapse to higher genus asymptotically AdS black holes. In a 
related work, Lemos [13] matched a flat FLRW metric to a radiating Vaidya 
exterior. Since their introduction by Lemos [14], the physics of toroidal and 
higher genus black holes has also been studied (see e.g. [5, 22, 15, 16]). 
Subsequently there has been considerable interest in toroidal and higher 
genus models of black holes, shells and horizons (see e.g. [10]), partly due 
to the existence of a "landscape" of vacua states in string theory with A 
positive, negative and zero (see e.g. [8]). 

In this paper we shall investigate generalisations of some of the above 
results by matching non-spherical inhomogeneous and anisotropic interiors 
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(mostly taken from the excellent book of Krasiriski [11]) to the non-spherical 
Kottler metric with A < 0, thus describing a wider class of examples of 
inhomogeneous and anisotropic gravitational collapse to toroidal and higher 
genus black holes. 

The Kottler metric with A < has initial and final curvature singular- 
ities and a time-like null infinity, J^. Any radial time-like geodesic has an 
endpoint at each singularity, crossing the future horizon of J? on its way 
to a maximum radius before falling back through the past horizon of J* to 
r = (see Figure 1). In a matching to a dust interior, such a geodesic is 
the outer boundary of the matter. The interior is part of a cosmological so- 
lution which in turn typically has initial and final singularities. As we shall 
see below, for the matching to FLRW, the initial and final singularities of 
the interior are matched to the initial and final singularities of the exterior, 
but this isn't necessarily true for all the matchings we find. However we may 
arrange that the initial singularities coincide and that the final singularity 
of the interior is behind the event horizon. 

In what follows we use units such that the speed of light c and Newton's 
gravitational constant G satisfy c = 8ttG = 1. 



2 Matching Theory in Brief 

Let (M^^g^) be spacetimes with non-null boundaries . Matching them 
requires an identification of the boundaries, i.e. a pair of embeddings $± : 
a — ► with 3>±(<t) = cr 1 * 1 , where a is an abstract copy of either boundary. 
Let be a coordinate system on a. Tangent vectors to cr^ are obtained by 
f i a = andHigherGenus-g^, and we shall usually work with orthonormal 

combinations ef a of the f^ a - There are also unique (up to orientation) 
unit normal vectors n± to the boundaries. We choose them so that if n" 
points into M + then n°L points out of M~ or viceversa. The first and second 
fundamental forms are simply qf- = ef a ef^ 'g a g\ . , Hf 1 - = —n^ef^V^ef 01 . 
It is well know that the matching conditions (in the absence of shells) require 
the equality of the first and second fundamental forms on a^, i.e. 

(1) q± = qT., H±=Hr.. 

These will be referred to as first and second matching conditions, respec- 
tively. If symmetries are present, one chooses the e^ a to reflect the symme- 
try and then iJ-- simplifies, as we shall see. 

It is a familiar fact that, in order to match a perfect fluid spacetime to a vac- 
uum solution, the fluid's pressure has to vanish at the matching boundary. 
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When the matching conditions are satisfied, this necessarily follows from 
the Hamiltonian constraint (on the time- like boundary). It is possible, but 
perhaps rather artificial, to match a perfect fluid interior with cosmological 
constant to a vacuum exterior with a different cosmological constant at a 
surface with non-zero pressure: the jump in the cosmological constant is 
chosen to match the jump in the pressure. However, we shall mostly be con- 
cerned with dust interiors with cosmological constant and vacuum exteriors 
with the same cosmological constant. The time-like boundary of the inte- 
rior must be ruled by matter world lines, which for dust are geodesies, and 
the matching conditions require that the boundary is also ruled by time-like 
geodesies of the exterior. This is part of the condition (1) matching the 
second fundamental forms: if we take the tangent to the fluid world lines as 
e^ Q then the geodesic condition implies the vanishing of H^. If the interior 
is not dust, as for example with the inhomogeneous p ^ perfect fluids of 
section 3.3.3, the matter world-lines ruling the boundary are accelerating 
and now the accelerations must match. 



3 Spacetimes and matching 



3.1 The Kottler solution as the exterior 



To begin with, let us recall the Kottler solution which we shall always take 
as the exterior in what follows. The metric is 

(2) ds\ = -VdT 2 + V~ l dr 2 + r 2 {d6 2 + Y J 2 {6)d^ 2 ) 

with 

2m A , 

3 V = b -r 2 

r 3 

and = 6, sin#, sinh#, according to b = 0,1,-1, respectively. With 

A < and m > 0, V has a unique positive zero and this solution describes 
a black hole with planar, spherical or hyperbolic symmetry respectively, on 
an asymptotically AdS background. For b = or —1 it is possible to make 
identifications of the 2-metric of constant T and r to obtain toroidal and 
higher genus black holes. Notice that will have the same topology (times 
R). 

Radial time-like geodesies of the metric (2) are characterised by the system 

VT = E 
r 2 = E 2 - V 

where E is a constant which is strictly positive for a future-directed time-like 
geodesic which is ever outside the black hole. Therefore, with A < 0, every 
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such geodesic has a maximum value of r, and so has two end-points, past 
and future, at r = (cf. Figure 1). For the dust interiors which follow, these 
geodesies are the candidate outer boundaries. 

r = (future singularity) 



horizon 




geodesic 



(past singularity) 



Figure 1. Penrose diagram for the Kottler solution with 
A < 0, showing a typical radial time-like geodesic at which 
the solution can be matched to an FLRW interior. The hori- 
zon is the event horizon for the J? at the right. 



(a) Toroidal Kottler: 

In the case b = 0, we introduce Cartesian coordinates (y, z) in place of 
(Q,4>)- To perform the matching preserving the symmetry we shall take a 
3-surface a + in the exterior generated by the orthonormal tangent vectors 

= Tdx + rd r , e\ = r~ 1 d y , = r~ 1 d z 

with ef tangent to geodesies, and parametrized by 

<S> + (\,y,S) = (T(\),r(\),y,S). 

An over dot will mean differentiation with respect to A, which is proper time 
along the geodesies ruling the boundary. Thus in particular 

(4) Vf 2 - y-V 2 = 1, 
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which, together with the constancy of TV, is equivalent to the geodesic 
equation. The induced metric on a + is 

(5) ds\\ a = -dX 2 + r 2 (dy 2 + dz 2 ). 

In all cases below, the normal from the exterior will be taken to be 

(6) n+ = yd T + Vfd r 

and the only non-zero terms of the second fundamental form are 

( 7 ) H 22 = -^33 = 



(b) Higher Genus Kottler: 

With b = —1, we can take the hyperbolic metric on the surfaces of con- 
stant (T, r) in the upper-half-plane form, so that the spacetime metric is 

r 2 

dsl = -Vdt 2 + V~ x dr 2 + -^{dy 2 + dz 2 ). 

yl 

The induced metric on the matching surface is 

2 

(8) ds 2 + L = -d\ 2 + r -^(dy 2 + dz 2 ). 

yl 

The normal to the matching surface can be written as in (6) and the non- 
zero components of the second fundamental in this case are again given by 
(7). 

In the next sections, we shall describe the different families of collapsing 
interior spacetimes to be considered, as well as their matching to the Kot- 
tler exterior. We shall divide our analysis into two parts, corresponding to 
spatially homogeneous and inhomogeneous collapsing spacetimes. 



3.2 Spatially homogeneous collapsing spacetimes as interiors 



3.2.1 FLRW spacetimes 

It is convenient to write the FLRW solutions as: 

(9) ds 2 _ = -dt 2 + R 2 (t){dp 2 + f 2 (p){d8 2 + g 2 (6)d<p 2 )), 

with the functions / and g as given in Table 1, for the different values of the 
curvature k. The cases with k = 1 all correspond to spherical symmetry. For 
k = or — 1, cases (a) have planar symmetry, (b) have spherical symmetry 
and (c) have hyperbolic symmetry. 
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k = 1 


k = 


k = -l 


/ = smp, cos p; 
g = sin 9, cos 9 


(a) / = 1 = g 

(b) / = p; g = sin 9, cos 9 


(a) / = e ±p ; g = 1 

(b) / = sinh p; g = sin 9, cos 9 

(c) / = cosh p; g = sinh 9, cosh 9 



Table 1. The / and g functions of FLRW metric (9) de- 
pending on the curvature k. 



As we noted in the introduction, the matching of a spherical form of the 
FLRW metric to a spherical exterior is familiar. The possible matchings 
of the non-spherical forms of FLRW metrics as interior to a non-spherical 
exterior Kottler solution can be summarized as follows: 

(i) k = and 6 = 0: Cannot match. This fact was given in [20] and 
[13]. In the latter reference, this FLRW was matched to a 6 = Vaidya 
exterior instead. 

(ii) k = -1 and b = -1: Matched in [20]. 
(hi) k = — 1 and 6 = 0: Can match. 

The third possibility was mentioned but not carried out in [20]. To see 
how it works, take the FLRW from the third column of Table 1 with / = e p 
and g = 1 and take a~ to be p = pq. The metric on the boundary is 

(10) ds 2 _\ a = -dt 2 + R 2 (t)e 2p0 (d9 2 + d^ 2 ), 
which must match to (5), so that (t,9,(f>) = (X,y,z) and 

(11) r = R(t)e P0 , 

where we introduce the notation == for an equality holding at the match- 
ing surface. The second fundamental form from the interior has nonzero 
components just 

H--H-- 1 

22 ~ 33 ~ R(ty 

which must match (7), so that 

(12) Vf = e p0 . 

The matching will be possible if the proper time dependence imposed on r 
by (11) is consistent with it defining a geodesic. Since (12) holds, this is 
equivalent to the normalisation condition (4), which, with the aid of (11) 
and (12), becomes 

(13) R 2 = 1 + - + -R 2 . 
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This can be recognised as the Friedmann equation for the interior, with a 
dust density p satisfying 

(14) -R 3 p = me- 3p0 . 

Thus (14) gives the mass m of the exterior in terms of the matter in the 
interior. Every solution of (13) that expands from an initial singularity has 
a maximum value of R and then recollapses, just as the time-like geodesies 
in the exterior do. By (11), the matching runs all the way from bang to 
crunch. Note that, even with the boundary identified as a torus, the interior 
is not spatially compact, as the range of p goes to — oo. 



3.2.2 Anisotropic spacetimes as interiors 

As an example of a homogeneous anisotropic interior, we consider the tilted 
dust LRS Bianchi V metric given by the Farnsworth solution [9, 21], which 
reduces to case (hi) of the previous section when the tilt vanishes. The 
metric is: 

(15) ds 2 _ = -dt 2 + f 2 dp 2 + g 2 e 2p (dy 2 + dz 2 ). 

Here / and g are functions of a single variable u = t + ap, a and a are 
constants and 

(16) f = a(ag u + g), 
while g satisfies 

lgg U u + g 2 - Ag 2 — ? = 0, 



which has the first integral 

(17) gg u = -g 6 + — + c 

O Ob 



with c another constant (compare (17) with (13)). The fluid flows along the 
vector field u = d/dt but the surfaces of constant t are not the surfaces of 
homogeneity: this is what it means for the fluid to be tilted. The energy- 
density is given by 

6(1 Ag - g uu ) 

V = i 

ag u + g 

which, with the aid of (17), simplifies to 

3ac 
1? 

The signs of the constants need to be chosen so that / and p are positive. 
As we shall see below, we need c positive, so also a > 0. Solutions of (17) 
expand from g = at say u = 0, to a maximum value g m ax an d then contract 
back to g = at say u = up. Near the beginning g ~ u 2 / 3 while near the end 



(18) p 
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g ~ (up — u) 2 / 3 . From (16), if we choose a > 0, then / ~ -u -1 / 3 near u = 0, 
with a positive constant of proportionality, but then / would be negative 
near u = uf, and so must vanish before u = uf- This will be a curvature 
singularity (since u will diverge). It must happen in the contracting phase 
(so that g u < 0) and will occur when g is the unique positive root of the 
cubic 

(19) 0(9):= (|_i_) / + J< +c , 

obtained by equating / to zero and using (17). It will be important, once 
we have matched this interior to the planar Kottler exterior, that this sin- 
gularity lies inside the event horizon. 

The Farnsworth solution generalises the k = — 1 FLRW solution which can 
be recovered for a = 1 and a = 0. We shall now show that this solution can 
be matched to Kottler with 6 = 0. The matching surface c~ is parametrized 
by 

$ {\y,z) = (X,p ,y,z). 
The first fundamental form in this case is 

ds 2 -\ a = -d\ 2 + g 2 e 2p0 (df + dz 2 ) 

so the first matching condition is 

(20) r=ge p0 . 
The unit normal to the matching surface is 

n = jdp 

and the second fundamental form is found to have non-zero components 

1 ( 9u 

so that the second matching condition is 



H22 = H 33 = - I a + 1 



VT a 1 / g 

- a h 1 



r f V 9 
which, with (20), implies 

yr = — , 

a 

a constant, as required. For (4) to hold, we need 

2 1 A 2 2me- 3p0 

9u = — + T9 H 

a 2 3 g 

which is the remaining field equation (17) if we make the identification 

c = 2me- 3p0 , 
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so that indeed c > 0. 

This completes the matching, which will hold from the initial singularity 
until / vanishes. Then it fails and the exterior can no longer be the Kottler 
solution. However, the vanishing of / can happen inside the event horizon, 
when the exterior solution will be Kottler up to and beyond the event horizon 
(see Figure 2). 

The vanishing of / happens when g is a root of (19). This will be inside 
the horizon if G(gn) < 0, where gn is the value of g corresponding to the 

1 /3 

location of the horizon, r = rjj = (— ^p) • The condition is 



which is therefore a condition on the parameters of the interior solution. It 
is easy to check that this condition delays the final singularity in the interior 
until after the formation of a marginally-outer-trapped surface (which occurs 
when g u = —a" 1 , or, using (17), when g = gjj). 



(21) 




r = (Kottler future singularity) 



i 



+ 



/ = (Farnsworth future singularity) 





Cauchy horizon 



horizon 



apparent horizon 



p = — oo 



matching surface 




past singularity 



Figure 2. Penrose diagram for the matching of Farnsworth 
to Kottler with A < 0. The matchings for Szekeres interiors 
are qualitatively similar. 
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3.3 Inhomogeneous collapsing spacetimes as interiors 



3.3.1 A higher genus analogue of the Farnsworth solution 

A dust solution (in the Szekeres family) which is a higher genus analogue 
of the Farnsworth solution, but not itself homogeneous, may be constructed 
as follows. The metric (compare (15)) is 

(22) ds 2 = -dt 2 + fdp 2 + g 2 cosh 2 p ^y 2 + dz2 ) 

where / and g are functions of a variable u = t + a log cosh p for constant a. 
Then the Einstein equations are satisfied for dust with 4- velocity d/dt and 
cosmological constant A if 

/ = ag u + g, 

and g satisfies 

2 A 3 

gg u = 3-5 +g + c , 

in terms of another (positive) constant c. The dust density is 

3c 
J? 



This solution is not homogeneous on the surfaces of constant u. It can be 
matched to the hyperbolic Kottler solution at p = po if 

r = g cosh po 
VT = sinh p 
when the normalisation condition (4) reduces to 

c = 2m cosh -3 po- 

Just as in the previous case, with a positive a there exist solutions which 
expand from a singularity at which g vanishes and / diverges, pass through 
a maximum g, and recontract. In the contracting phase, / will vanish before 
g and this will mark the final singularity. For this to lie inside the future 
horizon there is a condition analogous to (21), which this time is 

a < gn coth pq. 

The interpretation is as before: this condition ensures that a marginally- 
outer-trapped surface forms on the boundary before the interior becomes 
singular. A difference from the previous case is that po appears explicitly 
in the condition, and now, given the parameters (c, a, A) determining the 
interior, there will always be a (small enough) po at which to make the 
matching. 
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3.3.2 The Szekeres solutions 

We next consider the Szekeres form of the metric 

ds 2 = -dt 2 + e 2a dp 2 + e 2fS {dy 2 + dz 2 ) 

where a and (3 are two C 2 arbitrary functions of t, p, y, z. This is compatible 
with a perfect fluid source with 4- velocity u l = <5q and pressure p = p{t). 
We only consider the case where 8(3/ dp ^ 0, which is the Szekeres class that 
generalises the Lemaitre-Tolman spherically-symmetric dust solution. This 
is because we need the 2-surface of constant (t, p) to have non-zero second 
fundamental form in the 3-surface of constant t. 

The Einstein field equations for a perfect fluid imply 

(23) e a = h(p)e^d p e 13 

(24) e- v = A(p)(y 2 + z 2 ) + 2B 1 (p)y + 2B 2 (p)z + C(p) 

where A,B\,B2,C and h are arbitrary C 2 functions of p. The 2-surfaces of 
constant (t, p) have a metric which is conformal to 

e 2u (dy 2 + dz 2 ), 

with a conformal factor independent of y and z. From the form (24) of u, 
this 2-surface has constant positive, zero or negative curvature according as 
AC — B 2 — B\ is positive, zero or negative and we can conveniently call these 
the spherical, flat or hyperbolic cases respectively. 

To obtain a dust source with cosmological constant, we set p(t) = —A 
and then from [3] or [11] we find 

(j) 1 = -k(p) + — + -A(j) 2 



3 



3/M 



d p (2Me 
e 2 Pd p (eP) 



/' 
with 

(25) AC -B 2 -B 2 = \(h- 2 {p) + k{p)). 

Note the scaling freedom: 

(A,Bi,C) = e<^(A,Bi,C) 

v = v — £ 

h = e~ c /i 

k = e^k 
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and the freedom to perform affine transformations and inversions in y and 
z. 

The spherical case has already been matched to the spherical Kottler 
solution [12], so here we shall perform the matching of the flat (including 
toroidal) and hyperbolic (including higher genus) Szekeres solutions to the 
corresponding Kottler solutions. We shall match across p = po, and use 
affine transformations in y and z to simplify A, Bi and C at po. 

It is worth noting that the matchings we obtain for the non-spherical 
Szekeres solutions specialise to give matchings of the non-spherical Lemaitre- 
Tolman solutions (for these solutions see [11]). It is also worth noting that 
the Farnsworth solution of section 3.2.2 and its analogue of section 3.3.1 
fall into the Szekeres class of solutions, respectively planar and hyperbolic. 
Consequently one should expect that, as in those examples, there will be 
restrictions on the range of time over which the matchings that we find will 
hold, arising as conditions on the interior. 



(a) Toroidal Szekeres: 



Without loss of generality, we use the freedoms noted above to choose 
A = Bi = 0, C = 1 at po- The first fundamental form on the matching 
surface can be written as 

(26) ds 2 _\ a = -d\ 2 + 4> 2 (df + dz 2 ) 
so that the first matching condition gives 

(27) 0(A,p)=r. 
The normal to the matching surface is 

n~ = e~ a d p , 

and the second fundamental form at a has non-zero components 

1 

so that the second matching condition is 

(29) Vt = -. 

h 

From (4) we obtain the matching of the mass: 

M(p ) =m. 

This completes the matching. With M(p) and M tP (p) both positive (as is 
physically natural), the density is singular where either of 4> or /3 p vanishes. 
Vanishing (j) is matched, via (27), to the singularity in the exterior, but 
vanishing >p entails vanishing of e a via (23). This is like the problem of 



(28) H 22 - H 33 
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vanishing / in the previous two examples but, because of the extra functional 
freedom in the general Szekeres solution, we cannot write down an explicit 
condition like (21) to force this to happen behind the horizon. Rather we 
simply insist, as a condition on the interior, that /3 tP be positive for values 
of (j> giving values of r outside the horizon. 

(b) Higher genus Szekeres: 

By suitable transformations of (y, z) and scaling, we may assume without 
loss of generality that A = B2 = C = 0, B\ = 1/2 at po- The induced metric 
on the boundary can be written as 

ds 2 _L = -d\ 2 + tr(dy 2 + dz 2 ) 

and the first matching condition is 

(30) = r. 
The unit normal to a is 

rT = e~ a d p 
so that the second matching condition is 

VT = i 
h 

The normalisation condition (4) reduces to 

M(p ) = m 
which is continuity of the mass again. 

The comments at the end of the previous section regarding the range over 
which the matching holds apply in this case also. 

3.3.3 Barnes solutions 

The Barnes solutions [2] can be written in the form [11] 

ds 2 = -D 2 dt 2 + W~ 2 {dx 2 + dy 2 + dz 2 ) 

where W = W(t, x, y, z) is determined, in a way which varies from case to 
case, by another function w = w(t,p), with p = p(x,y,z), satisfying 

(31) ^ = f{p) 

w 

and 

W t 

D = F — - 
W 

with F = F(t) an arbitrary function. These solutions, which include the 
well-known McVittie class, may be characterised as the general solution of 
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the Einstein equations with the perfect-fluid form of stress-tensor, with a 
twist-free, shear-free but expanding 4-velocity. In general they don't have a 
barotropic equation of state (see the account in [11]). A solution is deter- 
mined by a choice of F(t) and f(p), with data w(t, po) and w p (t, po) for (31) 
at some value po of p. 

We require as before that the boundary be ruled by matter flow-lines, 
and that the pressure vanishes at the boundary, but since these are not dust 
solutions, we don't require that the matter flow-lines ruling the boundary 
be geodesies. 

Some examples of matching in the spherical case have been done already 
(see e.g [6] and other references in [11]) so we shall again consider only the 
plane and hyperbolic cases. 

(a) Toroidal Barnes: 

These solutions have 

W = w(t, p), with p = x 

and 

(32) p + A = JL+2fw 2 -3w 2 p 

(33) p-A = -j^ + w 2 p -2-^ -2^-(ww >tp -w tt w >p ). 

We take the matching surface from the exterior to be given in terms of 
proper time A by (T(A), r(A)) and from the interior by (t(X),p = po). The 
first matching conditions are 

(34) VT 2 - V- l r 2 ± D 2 i 2 ^ 1 



and 
(35) 

From (35) we obtain 
(36) 



1 

w 



1 

The normal to the matching surface is 



n = wd p 



and the second fundamental form gives 

H 22 = H 33 = -DwJ. 
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The second matching conditions are 

(37) -w p = — 

(38) ~ W ~D " Tr-rf- ^VV r f 3 + ^V^Vrfr 2 

where the extra complexity of (38) reflects the fact that the boundary is not 
ruled by geodesies; the right-hand-side in (38) is (minus) the inner product of 
the acceleration in the exterior of the boundary curve with the unit normal. 

The normalisation condition (4), with (36) and (37), leads to 

2 

( 39 ) v = " ivm 

and then preservation of this in time implies, by (33), 

(40) p = 

which also results from (38), given the rest of the matching. Thus, as ex- 
pected, the pressure must vanish at the matching surface. 

One may ask, given a planar Barnes solution, where should the match- 
ing surface be? We need a (necessarily comoving) surface of constant p, say 
p = po, at which p = 0. Then the quantity Q, defined by 

n 1 ( 1 2 A 

is necessarily constant on p = po, by virtue of the field equations. We match 
to an exterior planar Kottler solution according to (34) and (35) and then 
we satisfy (39) provided the parameter m of the Kottler exterior equals Q. 

Not every Barnes solution will contain a constant-p surface at which p 
vanishes. However it is possible to construct examples by using the match- 
ing surface as the surface at which data for (31) is given. That is, we choose 
F(t) and f{p), and, then given a choice of boundary p = po, we choose 
w(t,po) (positive) and w p (t,po) so that the density p(t,po) from (33) is pos- 
itive, the pressure p(t, po) from (33) is zero, and the quantity Q(po), which is 
now necessarily constant in time, is positive (we also want D positive and, to 
be matching outside the black hole, w tP negative). Integration of (31) then 
gives a Barnes solution, at least in a neighbourhood of the boundary, which 
has a boundary with the desired properties. In the exterior, the boundary 
follows a curve given by (35) together with the integration of (37), using 
i = D-K 



(b) Higher genus Barnes: 
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In this case 

W = yw{t,p) with p = x/y 
and therefore the metric can be written as 

ds 2 = -D 2 dt 2 + -^-^(y 2 dp 2 + 2pydpdy + (1 + p 2 )dy 2 + dz 2 ) 
y l w l 

where we have changed variables to eliminate x, together with 

(41) n + A = + 2(p 2 + l)/u; 3 + 6pw, p - 3(p 2 + l)w 2 p - 3w 2 

3 „ . , 9 . s 2 . 2 o 



P-A = -7; ^ - 2pura; iP + (p^ + 1 W p + w z -2- 



(42) +2^ (p - (p 2 + 1)^) («;«, )tp - 

We take the matching surface from the exterior to be given by (T(A),r(A)) 
and from the interior (i(A),p = po). The first matching conditions are 

(43) Vf 2 - V- X r 2 = D 2 i 2 = 1 
and 



(44) r 2 * l±P 

The normal to the matching surface is 



2 



n = w 



[7-. — op> w py « 

V l + P 



and the second fundamental form is 

:D 
~D 



H n = -w^l + p 2 ^- 



H 22 = H 33 = i^—a ~ Vl + P 2 W,p 



so that the second matching conditions are 

(45) p-(l + p 2 ) ^P * V f 
w 



~>D n a " . .. ■ 1 T , T , -q 3, 



(46) -wy/l + P 2 ~jJ- = fr-rf--VV r f 3 + -V~ 1 V r fr 2 . 
It follows from (43) and (45) that 

v=(p-(i+p 2 )^y-(i+p 2 )-±- 2 , 

and then it is straightforward to confirm that 

p = 0. 

One argues as before for the existence of higher-genus Barnes solutions with 
a constant-p surface at which p vanishes: choose the free functions F(t), 
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f(p) and the data w(t, po), w iP (t, po) so that p(t, po) = 0; propagate inwards 
with (31) to obtain the interior; the quantity 

( p 2 + 1) l/2/ A(p 2 + 1) (1+p 2 ) 2 2 \ 

Q = 2u;3 " 3 + ~ (P W ~ (P + 1 ^ J 

is then necessarily constant in time at po and under the matching equates to 
m in the exterior. The data and free functions need to satisfy various open 
conditions to ensure that m, V, T and so on are positive but this is easily 
done. 
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